7  Power  Sou 


213  (2012)  16-30 


Contents  lists  available  at  SciVerse  ScienceDirect 

Journal  of  Power  Sources 


ELSEVIER 


epage:  www.e 


■n/locate/jpowsou 


Alkaline  membrane  fuel  cell  (AMFC)  modeling  and  experimental  validation 

E.M.  Sommer3*3,  L.S.  Martins'3,  J.V.C.  Vargas3  *,  J.E.F.C.  Gardolinski c,  J.C.  Ordonez13,  C.E.B.  Marino3 

1  Departamento  de  Engenharia  Mecanica,  Universidade  Federal  do  Parana,  C.  P.  19011,  Curitiba,  Parana  81531-990,  Brazil 
b  Department  of  Mechanical  Engineering  and  Center  for  Advanced  Power  Systems,  Florida  State  University,  Tallahassee,  FL  32310-6046,  USA 
cPANalytical  B.V.,  Lelyweg  1,  7602EA,  Almelo,  The  Netherlands 


ARTICLE 


N  F  O 


A 


S  T  R  A  C  T 


Article  history: 

Received  3  March  2012 
Received  in  revised  form 
22  March  2012 
Accepted  23  March  2012 
Available  online  11  April  2012 


Keywords: 

Alkaline  membrane  fuel  cell 
Potassium  hydroxide 
Dynamic  mathematical  modeling 
Experimental  validation 


This  paper  aims  to  produce  a  dynamic  model  that  is  computationally  fast  to  predict  the  response  of  the 
single  AMFC  according  to  variations  of  physical  properties  of  the  materials,  and  operating  and  design 
parameters.  The  model  is  based  on  electrochemical  principles,  and  mass,  momentum,  energy  and  species 
conservation.  It  also  takes  into  account  pressure  drop  in  the  gas  channels  and  the  temperature  gradient 
with  respect  to  space  in  the  flow  direction.  The  simulation  results  comprise  temperature  distribution,  net 
power  and  polarization  curves,  which  were  experimentally  validated  by  direct  comparison  to  voltage  and 
current  measurements  performed  in  a  cellulose-based  AMFC  prototype  for  different  electrolyte  (KOH) 
solution  concentrations  (y),  showing  good  quantitative  and  qualitative  agreement.  It  is  concluded  that 
the  startup  transient  is  short  and  that  there  are  optimal  values  of  y  ( ~40  wt.  %  )  which  lead  to  maximum 
power,  that  are  herein  shown  experimentally  for  the  first  time.  In  the  process,  the  model  was  used  to 
formulate  empirical  correlations  for  the  exchange  current  density  (io)  in  the  electrodes  with  respect  to 
the  electrolyte  concentration  for  future  fuel  cell  development.  Therefore,  the  adjusted  and  validated 
model  is  expected  to  be  a  useful  tool  for  AMFC  control,  design  and  optimization  purposes. 

©  2012  Elsevier  B.V.  Oven  access  under  the  Elsevier  O  A  license. 


1.  Introduction 

Fuel  cell  technology  is  currently  at  an  advanced  stage.  However, 
there  are  still  challenging  hurdles  to  be  overcome  to  allow  for 
extensive  utilization,  such  as  improvements  on  hydrogen  produc¬ 
tion,  storage  and  distribution,  catalysts  cost  reduction  and 
manufacturing  [1,2], 

Among  the  existing  fuel  cell  types  [2],  Alkaline  Fuel  Cells  (AFC) 
have  interesting  features  as  compared  to  polymer  electrolyte  fuel 
cells  (PEMFC),  such  as  higher  current  density,  lower  cost  electrolyte 
(KOH  aqueous  solution)  and  mainly,  the  possibility  of  using  non¬ 
noble  catalysts  (e.g.,  nickel,  silver),  since  faster  kinetics  of  the 
reactions  is  observed  in  alkaline  media  than  in  acid  media  [3,4], 
Additionally,  AFC  can  work  at  higher  temperatures  (100-120  °C) 
than  PEMFC,  which  cannot  operate  above  90  °C  due  to  the  need 
to  hydrate  the  naflon  membrane  [4],  However,  PEMFC  are  more 
suitable  for  vehicular  applications  mainly  due  to  the  use  of  solid 
electrolytes  that  eliminate  the  risk  of  leakage  [2], 

The  electrolyte  in  the  AFC  consists  of  an  aqueous  alkaline 
solution,  such  as  potassium  hydroxide  (KOH).  Since  aqueous  alka¬ 
line  solutions  absorb  carbon  dioxide  (CO2),  an  undesirable 
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"poisoning"  effect  is  observed  through  the  conversion  of  KOH  to 
potassium  carbonate  (K2CO3),  in  a  reaction  that  competes  with  the 
expected  redox  reaction  between  hydrogen  and  oxygen  at  the 
anode  and  cathode  [4],  This  reaction  has  the  effect  of  reducing  the 
number  of  hydroxyl  ions  available  for  reaction  at  the  electrodes, 
also  reducing  the  ionic  conductivity  of  the  electrolyte  solution,  and 
could  block  the  pores  of  the  gas  diffusion  layer  (GDL)  by  the 
precipitation  of  K2CO3  salt,  but  not  causing  any  degradation  of 
the  electrodes  [5],  As  a  result,  the  AFC  performance  is  greatly 
reduced  in  air  breathing  systems,  i.e.,  when  the  oxidant  is  air  which 
contains  CO2. 

Because  of  the  CO2  poisoning  effect,  alkaline  fuel  cells  typically 
operate  on  pure-oxygen,  or  at  least  purified  air,  and  must  have 
built-in  mechanisms  to  clean  out  as  much  of  the  carbon  dioxide  as 
is  possible.  Because  the  generation  and  storage  requirements  of 
oxygen  make  pure-oxygen  AFCs  expensive,  the  active  development 
of  the  technology  has  not  attracted  large  industrial  interest. 
Although  the  harmful  effects  of  carbon  dioxide  poisoning  can  be 
partly  reduced  by  several  different  strategies  [4],  a  permanent 
solution  would  increase  the  possibility  of  AFC  commercialization. 

In  order  to  cope  with  that  reality,  two  types  of  AFCs  exist:  static 
electrolyte  and  flowing  electrolyte.  Static,  or  immobilized,  elec¬ 
trolyte  cells  are  of  the  type  used  in  the  Apollo  space  craft  and  the 
space  shuttle,  which  typically  use  an  asbestos  separator  saturated 
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Nomenclature 

Ac  total  gas  channel  cross-section  area,  m2 
As  fuel  cell  cross-section  area,  m2 

Aw  wall  heat  transfer  area,  m2 

AFC  alkaline  fuel  cell 

AMFC  alkaline  membrane  fuel  cell 

B  dimensionless  constant 

Ba  bias  limit  of  quantity  a 

C  constant 

Q  concentration  of  oxidized  and  reduced  states 

(i  =  0,  R),  mol  m-3 

CFD  computational  fluid  dynamics 

CVi  control  volume  i 

c  specific  heat,  kmol  kg-1  K-1 

cp  specific  heat  at  constant  pressure,  kj  kg  1  K-1 

D  Knudsen  diffusion  coefficient,  m2  s-1 

Dh  gas  channel  hydraulic  diameter,  m 

/  friction  factor 

F  Faraday  constant,  96,500  C  eq-1 

h  heat  transfer  coeficient,  W  m-2  k-1 

Hj(Ti)  molar  enthalpy  of  formation  at  temperature  Ti  of 
reactants  and  products,  kj  kmol-1  of  compound  i 
Hj(Oj)  dimensionless  molar  enthalpy  of  formation 

at  temperature  0;  of  reactants  and  products  of 
compound  i 

ippe  inverse  problem  of  parameter  estimation 

!*o,a>  io.c  anode  and  cathode  exchange  current  densities,  A  m-2 

iLim,a,  ium,c  anode  and  cathode  limiting  current  densities,  A  m-2 

/  current,  A 

j  mass  flux,  kg  m2  s-1 

k  thermal  conductivity,  W  m-1  K-1 

k°  standard  rate  constant,  ms-1 

K  permeability,  m2 

L  control  volume  length,  m 

Lc,  Lt  gas  channels  internal  dimensions,  m 

Lx,  fy,  Lz  fuel  cell  length,  width  and  height,  respectively,  m 

m  mass,  kg 

m  mass  flow  rate,  kg  s-1 

M  molecular  weight,  kg  kmol-1 

n  equivalent  electron  per  mole  of  reactant,  eq  mol-1 

h  molar  flow  rate,  mol  s-1 

N  dimensionless  global  wall  heat  transfer  coefficient 

nc  number  of  parallel  ducts  in  gas  channel 
p  pressure,  N  m-2 

P  dimensionless  pressure 

Pa  precision  limit  of  quantity  a 

ps  perimeter  of  cross-section,  m 

PEMFC  proton  exchange  membrane  fuel  cell 
Pr  Prandtl  number,  pcp/k 

q  tortuosity 

Q  reaction  quotient 

Q.  heat  transfer  rate,  W 

r  pore  radius,  m 

R  ideal  gas  constant,  kj  kg-1  K-1 

R  universal  gas  constant,  8.314  kj  kmol-1  K-1 

ReDh  Reynolds  number  based  on  Dh,  uDhp/p 
t  time,  s 

T  temperature,  K 

u  mean  velocity,  ms-1 

U  global  wall  heat  transfer  coefficient,  W  m-2  K-1 

lfa  uncertainty  of  quantity  a 


V  electrical  potential,  V 

VT  total  volume,  m3 

y  electrolyte  solution  mass  fraction,  wt.  % 

W  electrical  work,  J 

W  electrical  power,  W 

W  dimensionless  fuel  cell  electrical  power 

IV net  dimensionless  fuel  cell  net  power 

Wp  dimensionless  required  pumping  power 
[  •]  molar  concentration  of  a  substance  (mol  l-1) 

Greek  Symbols 

a  charge  transfer  coefficient 

fi  electrical  resistance,  Q 

y  ratio  of  specific  heats 

5  gas  channel  aspect  ratio 

AG  molar  Gibbs  free  energy  change,  kj  kmol-1H2 

AH  molar  enthalpy  change,  kj  kmol-1H2 

A P  dimensionless  pressure  drop  in  gas  channels 

AS  molar  entropy  change,  kj  kmol-1 

AT  temperature  change,  I< 

£  stoichiometric  ratio 

rja,  r]c  anode  and  cathode  charge  transfer  overpotentials,  V 
tjd ,a,  Vd,c  anode  and  cathode  mass  diffusion  overpotentials,  V 
rjohm  fuel  cell  total  ohmic  potential  loss,  V 

$  dimensionless  temperature 

H  viscosity,  kg  m-1  s-1 

Vi  reaction  coefficient 

if  dimensionless  length 

p  density,  kg  m-3 

a  electrical  conductivity,  fi-1  m-1 

%  dimensionless  time 

(p  porosity 

i/'  dimensionless  mass  flow  rate 

Subscripts 

a  anode 

(aq)  aqueous  solution 

c  channel  or  cathode 

cel  solid  part  of  electrolyte 

dif  diffusion 

e  reversible 

f  fuel 

(g)  gaseous  phase 

H2  hydrogen 

H20  water 

i  substance  or  location  in  AMFC 

i  irreversible 

i,a  irreversible  at  the  anode 

i,c  irreversible  at  the  cathode 

j  location  in  AMFC 

KOH  potassium  hydroxide 

(1)  liquid  phase 

max  maximum 

memb  electrolyte  membrane  (solid  and  liquid  parts) 

ohm  ohmic 

opt  optimal 

ox  oxidant 

OH-  hydroxide  ion 

02  oxygen 

ref  reference  level 

s,a  anode  solid  part 

s,c  cathode  solid  part 

s  solid  phase 
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electrolyte  solution 

wall 

wetted 

initial  condition 
control  volumes,  Fig.  2 
interaction  between  CV1  and  CV2 
interaction  between  CV2  and  CV3 
interaction  between  CV3  and  CV4 
interaction  between  CV4  and  CV5 


56  interaction  between  CV5  and  CV6 

67  interaction  between  CV6  and  CV7 

°°  ambient 

Superscript 

°  standard  conditions  [gases  at  1  atm,  25  °C,  species  ir 

solution  at  1  M,  where  M  is  the  molarity  =  (moles 
solute)/(liters  solution)] 

~  dimensionless  variable 


in  potassium  hydroxide,  and  more  recently,  the  use  of  solid  anionic 
membranes  which  are  meant  to  eliminate  the  carbon  dioxide 
problem  [6],  such  as  the  ammonium-based  Alkaline  Anion- 
Exchange  Membranes  (AAEMs)  [7],  In  the  flowing  electrolyte 
type,  a  sufficiently  open  matrix  is  used  to  allow  for  liquid  electro¬ 
lyte  circulation  to  reduce  poisoning  effects  [8],  and  the  old  elec¬ 
trolyte  may  be  exchanged  for  a  fresh  one,  as  the  AFC  becomes 
poisoned. 

Considering  static  electrolyte  cells,  the  use  of  asbestos  or 
ammonium-based  membranes  for  commercial  applications  also 
bring  additional  restrictions,  based  on  ammonia  potential 
hazardous  effects,  such  as  acute  toxicity  [9]  and  Immediately 
Dangerous  to  Life  or  Health  (IDLH)  levels  [10],  and  on  asbestos 
being  a  known  carcinogenic  agent  [11],  Therefore,  recently 
a  cellulose-based  alkaline  membrane  fuel  cell  was  proposed  [12], 
A  general  view  of  the  utilized  experimental  rig  and  of  the 
cellulose-based  membrane  AMFC  prototype  is  shown  in  Fig.  la 
and  b,  respectively.  Although  promising  experimental  results 
were  obtained  in  the  laboratory,  design  and  operating  parameters 
need  to  be  optimized  for  maximum  performance,  so  that  a  proper 
assessment  can  be  made  for  possible  future  commercialization. 
For  that,  mathematical  modeling  and  computational  simulation  is 
recommended  in  order  to  predict  how  fuel  cell  performance  is 
affected  by  the  variation  of  operating  and  geometric  parameters 
[13],  It  is  also  important  that  the  mathematical  model  accounts 
for  temperature  spatial  gradients  and  gas  channels  pressure 
drops  [13,14]  which  significantly  affect  fuel  cell  performance. 
Such  desired  model  features  should  not  lead  to  the  understanding 
that  the  more  complex  the  better  the  model  will  be.  On  the 
contrary,  when  the  model  is  considered  for  fuel  cell  structure 
optimization,  a  large  number  of  different  configurations  tests  are 
needed,  then  it  is  desirable  for  the  model  to  be  simple,  compu¬ 
tationally  fast  and  still  able  to  capture  reliably  the  fuel  cell 
behavior  [1,14], 

Proton  exchange  membrane  fuel  cells  (PEMFC)  and  solid  oxide 
fuel  cells  (SOFCs)  are  considered  the  most  promising  candidates  for 
commercial  exploitation,  and  much  attention  has  been  given 
toward  system  modeling,  including  the  fluid  flows  that  supply  the 
cells  with  fuel  and  oxidant,  with  the  use  of  computational  fluid 
dynamics  (CFD),  as  pointed  out  by  Young  [15]  and  Ma  et  al.  [16], 
that  recently  reviewed  the  status  of  three-dimensional  fuel  cell 
modeling.  On  the  other  hand,  there  are  few  previous  studies  con¬ 
cerning  AFC  modeling  and  simulation.  For  example,  Vargas  et  al.  [1  ] 
performed  the  AFC  internal  structure  thermodynamic  optimization 
introducing  a  dynamic  mathematical  model  that  accounted  for 
temperature  and  pressure  spatial  gradients;  Verhaert  et  al.  [17] 
presented  a  steady  state  AFC  thermodynamic  model  without  the 
consideration  of  possible  temperature  variation  within  the  cell 
using  a  control  volume  approach,  which  was  validated  by  experi¬ 
mental  data  for  future  fuel  cell  stack  development;  Duerr  et  al.  [18] 
used  an  isothermal  dynamic  electrochemical  model  to  predict  the 
effect  of  several  parameters  on  an  AFC  stack  response  during  load 
change  events;  Mohan  and  Shrestha  [19]  investigated  the  effect  of 
four  different  electrolyte  concentrations  and  three  different  anode 


flow  rates  on  AFC  performance  experimentally  and  with  an 
isothermal  steady  state  model,  concluding  that  performance 
increased  as  electrolyte  concentration  increased,  and  deduced 
statistically  from  modeling  and  statistical  analysis  optimal  elec¬ 
trolyte  concentration  and  anode  flow  rate  for  maximum  perfor¬ 
mance,  and  Weydahl  et  al.  [20]  introduced  an  isothermal,  one¬ 
dimensional  transient  model  of  a  porous  oxygen  electrode  in  an 
alkaline  medium  based  on  mass  balances  and  flooded-agglomerate 
theory  in  order  to  study  time  constants,  the  combined  effects  of  fast 
potential  propagation  and  slow  oxygen  diffusion  targeting  the 
dynamic  improvement  of  the  AFC  cathodes  without  sacrifice  of 
a  steady  state  response.  However,  no  studies  were  found  in  the 
open  literature  on  transient  mathematical  models  designed 
specifically  for  single  alkaline  membrane  fuel  cells  (AMFC)  that 
consider  all  of  their  internal  components. 

Though  very  useful  refined  information  can  be  extracted  from 
spatially  dependent  models,  the  two-  and  three-dimensional 


Fig.  1.  (a)  Fuel  cell  testing  worktop;  (b)  Alkaline  Membrane  Fuel  Cell  prototype. 
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models  are  not  suitable  for  the  optimization  of  flow  geometry, 
because  they  would  require  the  solving  of  partial  differential 
equations  for  flow  simulation  in  a  very  large  number  of  flow 
configurations.  Therefore,  in  this  paper  a  mathematical  model  is 
introduced  for  the  single  AMFC,  which  divides  the  fuel  cell  into 
several  control  volumes  that  correspond  to  the  most  representative 
parts  of  the  flow  system.  All  the  flow  phenomena  that  are  present 
are  taken  into  account.  The  result  is  a  model  with  unidirectional 
internal  flow  that  contains  additional  three-dimensional  features 
such  as  the  electrode  wetted  area,  heat  transfer  between  the  cell, 
fuel,  oxidant  and  the  surroundings,  and  pressure  drops  in  the  gas 
channels.  The  model  comprises  a  system  of  ordinary  differential 
equations,  the  solution  of  which  consists  of  the  temperatures  and 
pressures  of  each  control  volume,  and  the  polarization  and  net 
power  curves  for  the  whole  system.  The  model  is  simple  enough  to 
ensure  small  computational  time  requirements,  so  that  it  is 
possible  to  predict  the  response  of  the  single  AMFC  according  to 
variations  of  the  material  physical  properties,  and  operating  and 
design  parameters. 

The  mathematical  model  is  then  experimentally  validated 
using  5  (five)  different  sets  of  measured  data  obtained  for  the 
cellulose-based  AMFC  shown  in  Fig.  1  for  several  electrolyte 
concentrations.  First,  the  inverse  problem  of  parameters  estima¬ 
tion  methodology  is  used  to  elaborate  an  empirical  correlation  for 
the  electrodes  exchange  current  density  as  a  function  of  electro¬ 
lyte  concentration  using  3  (three)  data  sets,  therefore  adjusting 
the  mathematical  model.  After  that,  the  conditions  corresponding 
to  the  other  2  (two)  data  sets  are  computationally  simulated  with 
the  adjusted  model  to  verify  the  agreement  between  numerical 
and  experimental  results,  so  that  the  model  could  be  considered 
experimentally  validated  for  simulation,  control,  design  and  opti¬ 
mization  purposes. 


2.  Mathematical  model 

The  mathematical  model  is  based  on  an  AFC  model  introduced 
previously  [1],  with  the  inclusion  of  pressure  drop  in  the  gas 
channels,  heat  generated  by  overpotential  losses  and  the  cellulose- 
based  alkaline  membrane.  Although  the  single  AFC  mathematical 
model  has  been  presented  by  the  authors  before  [1],  considerable 
modifications  have  been  introduced  in  order  to  appropriately 
capture  the  AMFC  behavior.  Also,  the  nondimensionalization  of  the 
fuel  cell  geometry  has  been  improved  to  account  for  a  realistic  fuel 
cell  total  fixed  volume  engineering  design  constraint,  i.e.,  the  finite 
space  availability.  Hence,  in  order  to  provide  a  clear  understanding 
of  the  model  and  the  experimental  validation  procedure,  it  is 
instructive  to  summarize  its  main  features  together  with  the  new 
aspects. 

The  fuel  cell  is  divided  into  seven  control  volumes  (CV)  that 
interact  energetically  with  each  other  and  with  the  environment  as 
shown  schematically  in  Fig.  2.  The  seven  control  volumes  are:  fuel 
channel  (CV1),  the  anode  diffusion  layer  (CV2),  the  anode  reaction 
layer  (CV3),  the  cellulose-based  alkaline  membrane  (CV4),  the 
cathode  reaction  layer  (CV5),  the  cathode  diffusion  layer  (CV6)  and 
the  oxidant  channel  (CV7). 

The  model  consists  of  writing  mass,  energy  and  species 
conservation  equations  for  each  control  volume,  considering  the 
chemical  reactions  on  CV3  and  CV5.  All  flow  phenomena  are 
considered,  resulting  on  a  time  dependent  unidirectional  internal 
flow  model.  In  the  balance  of  energy,  the  heat  generated  by  the 
electrochemical  reactions  and  the  potential  losses  are  taken  into 
account.  The  heat  generation  mechanism  is  due  to  ohmic  resistance 
(throughout  the  fuel  cell),  as  well  as  by  activation  and  concentra¬ 
tion  overpotential  losses  in  CV3  and  CV5. 


The  solution  of  the  differential  and  algebraic  equations  yields 
temperature  and  pressure  profiles  for  each  control  volume,  and 
polarization  and  power  curves  of  the  system  as  well.  The  actual 
electrical  potential  and  power  are  obtained  by  subtracting  from  the 
reversible  potential  the  losses  due  to  surface  overpotentials  (poor 
electrocatalysis),  slow  diffusion  and  all  internal  ohmic  losses 
through  the  cell  (resistance  of  individual  cell  components, 
including  electrolyte  membrane,  bipolar  plates,  interconnects  and 
any  other  cell  components  through  which  electrons  flow).  These 
are  functions  of  cell  voltage  or  of  the  total  cell  current  (I),  which 
have  a  one-to-one  relationship  with  each  other.  Since  the  external 
load  determines  the  fuel  cell  operating  current,  in  this  treatment 
the  total  cell  current  is  considered  the  independent  variable.  For  the 
sake  of  generality,  and  for  control  and  real  time  operation  purposes, 
the  following  analysis  is  for  unsteady  state. 


2.1.  Fuel  cell  geometry 

Fig.  2  shows  all  the  external  geometric  parameters  of  the  fuel 
cell.  Lx  is  the  total  length  of  the  fuel  cell  (flow  direction),  Lj  the 
length  of  each  control  volume  of  the  AMFC,  such  that  1  <  j  >  7.  Ly 
and  Lz  are  the  height  and  the  thickness,  respectively. 

The  total  volume  of  the  fuel  cell,  Vj=  LxLy  Lz,  is  finite  and  fixed. 
This  is  a  realistic  design  constraint,  which  accounts  for  the  finite¬ 
ness  of  the  available  space  and  the  general  push  for  doing  the  most 
with  limited  resources  (e.g.,  space),  so  that  the  maximization  of 
performance  for  the  specified  volume  means  the  maximization  of 
performance  density.  The  fixed  length  scale  V^3  is  used  for  the 
purpose  of  nondimensionalizing  all  the  lengths  that  characterize 
the  fuel  cell  geometry,  as  follows: 


where  the  subscript  j  indicates  a  particular  dimension  of  the  fuel 
cell  geometry,  Fig.  2. 

The  wall  heat  transfer  area  of  one  control  volume  is  A™  =  ps  U 
(2  <  i  <  6)  and  A™  s  ps  L,-  +  tyIz  (i  =  1, 7;  assuming  that  Lt «  Lc  in 
Fig.  2),  where  ps  =  2  (Ly  +  Lz)  is  the  perimeter  of  the  fuel  cell  cross- 
section.  The  control  volumes  are  Vj  —  ty  Lz  Lj  (2  <j  <  6)  and  Vj  —  nc  Ic 
Lj  Lz  (j  -  1,  7),  where  nc  is  the  integer  part  of  Iy/(Lt  +  Lc),  i.e.,  the 
number  of  parallel  ducts  in  each  gas  channel  (fuel  and  oxidant). 


2.2.  Dimensionless  variables 


Dimensionless  variables  are  defined  based  on  system  geometry 
and  operating  parameters.  Temperature  and  pressure  are  refer¬ 
enced  to  ambient  conditions  (T„  and  p„),  —  Ti/T„  and 

Pj  =  Pi/Poo. 

Dimensionless  mass  flow  rate  (i/0,  global  wall  heat  transfer 
coefficient  (N),  heat  transfer  coeficient  (h),  thermal  conductivity  (k), 
area  (A),  specific  heat  (7),  and  time  (t)  are  defined  as  follows: 


=  mL  N  UwjVf *  h  hjVf  k  k^3 

rhref’  '  ^refcp,f  *  ^refcp,f '  *  ^refcp:f 

Aj  Cp  i  t 

=  77573:  Vi  =  tS  *  =  i- 


where  subscript  i  indicates  a  substance  or  a  location  in  the  fuel  cell, 
mre f  —  p„  Vj/(RfTaj  tref)  is  a  reference  mass  flow  rate,  and  Rf  the 
ideal  gas  constant  of  the  fuel,  and  tref  a  specified  reference  time. 
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2.3.  Mass  balance 


The  mass  and  energy  balances  for  CV1  state  that: 


The  hydrogen  mass  flow  rate  required  by  a  single  AMFC  for  the 
current  (I)  dictated  by  the  external  load  is 


mH2  = 


(3) 


Where  h  is  the  molar  flow  rate,  kmol  s-1;  M  the  molar  weight, 
kg  kmol-1;  n  the  equivalent  electron  per  mole  of  reactant, 
eq  mol-1 ;  /  the  total  current,  A,  and  F  the  Faraday  constant,  C  eq 
Accordingly,  the  oxygen  mass  flow  rate  needed  for  a  single  AMFC  is 
iho2  —  l/2fiH2Mo2.  For  simplicity  of  the  analysis,  it  is  hereby 
defined  that  ho2  —  2rho2/M02,  so  that  n02  =  nH2. 

2.4.  Energy  conservation 

2.4.1.  Fuel  and  Oxidant  Channels  (CV1  and  CV7) 

The  fuel  (hydrogen)  enters  the  system  through  the  CV1  channels 
and  the  oxidant  (oxygen)  through  the  CV7  channels,  as  shown 
in  Fig.  2. 


W  =  [^wl  +*f(8f-8i)+Q-n  +  Qiohm]  PfnJ^z  (4) 

where  6  is  the  dimensionless  temperature,  Q,  =  Q.i/threfCP]fT„  the 
dimensionless  heat  transfer  rate,  Q12  =  hiAs(l  -  q>2)(^2  -  #i). 
Qwi  —  NjAwi(l  -0j),  Q.j0hm  —  /^/(threfCp/Too),  /3  the  electrical 
resistance,  dig  the  initial  condition,  (ft  the  porosity,  As  =  LyLz/Vj'ri, 
yp f  =  si  V/h2  ■  and  ?  the  stoichiometric  ratio  which  is  assumed 
greater  than  1  on  both  sides. 

The  pressure  drop  in  the  channels  is  calculated  using  the  ideal 
gas  model  and  assuming  that  the  channels  are  straight  and  suffi¬ 
ciently  slender,  as  follows: 


where  i  =  1,  7  and  j  =  /,  ox,  respectively.  Here 

hi  =  (hjin  +  hi, out) /2  is  the  gas  dimensionless  mean  velocity  in 
the  channel,  defined  as  u  =  u/(RfTal)1/2,  f  is  the  friction  factor. 
According  to  mass  conservation,  the  flow  dimensionless  mean 
velocities  in  the  gas  channels  are  u i  =  C0\  (\ftf  -  i/,H2/2)/(AclPj), 
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u7  -  R0XO?7(i//0X  -  ip0J2)/(RfAc7Pox),  where  C  =  (RfT+)'/2mref/ 
(p„V?\  and  Aci  =  ncLcLilV^3,  i  =  1,  7,  is  the  dimensionless 
total  duct  cross-section  area  in  the  fuel  and  oxidant  channels, 
respectively. 

The  Reynolds  number  is  monitored  to  determine  flow  regimes. 
Appropriate  correlations  are  used  according  to  the  observed  flow 
regime  [21,22],  as  follows: 

Laminar  Regime  -  ReDh<2300  [21  ] 

fReDh.  =  24(l  -  1.35535}  +  1.94675?  -  1.70125?  +  0.95645f 
-  0.25375?) 

(6) 

hi^  =  7.541  (l  -  2.6105}  +  4.9705?  -  5.1 195?  +  2.7025? 

-  0.5045?)  (7) 

Turbulent  Regime  [22] 

/i  =  0.079ReD?/4  (2300  <  ReDhi  <  5  x  106)  (8) 


ftjPhj 

% 


(fi/gg>h,i  -  ?o3)Pn 

1  +  12.7(/}/2)1?2(Pr?/3  -  1) 


(2300  <  ReDhi  <  5 


106) 

(9) 


where  5f  =  Lc/Lh  Dhi  =  2LcL,/(Lt  +  Le)  for  i  =  1,7,  and 

ReDhi  =  U}Dh }/)}/ /*}. 

The  overall  fuel  cell  reaction  and  the  mass  balance  for  CV7  yield 
tiH2o  =  t&H2o,out  =  nH2o,in  =  ho2.  Therefore,  the  mass  and  energy 
balances  for  CV7  state  that: 


M  the  gas  molecular  weight,  q>  the  porosity,  q  the  tortuosity 
[1,24,25],  and  B  a  dimensionless  coefficient. 

The  output  pressure  of  the  diffusive  layer,  which  is  equal  to  the 
input  pressure  of  the  catalyst  layer,  is  obtained  using  Eq.  (11)  and 
the  ideal  gas  model,  as  follows: 


jiRiToo  kdi 

DiP  CO 


(12) 


where  J2  =  rhH2/A3wet  and  je  =  mQ2/A5wet,  and  A3,Wet  and  As,Wet  are 
the  wetted  areas  in  the  porous  catalyst  layers.  Note  also  that 
P2, in  =  Pf  and  P6,in  —  Pox-  The  average  pressures  in  CV2  and  CV6  are 
estimated  as  Pf  =  l/2(Pi  in  +  Pi  out),  i  =  2,6 

The  wetted  areas  in  the  porous  catalyst  layers  A3,wet  and  AsiWet 
are  estimated  by  considering  dual-porosity  sintered  metal  elec¬ 
trodes,  in  which  the  average  pore  diameter  is  considered  to  be 
roughly  the  same  order  as  the  square  root  of  the  porous  medium 
permeability  [1],  In  this  consideration  the  pores  are  approximated 
as  parallel  tubes,  i.e.,  Aj  wet  =  {4<pjLjKT1/2  +  (1  -  <pj)}As. 

Hydrogen  diffuses  through  CV1  and  CV3  and  oxygen  and  water 
vapor  diffuse  through  CV5,  CV6  and  CV7.  Therefore,  the  mass  and 
energy  balances  for  CV2  state  that 


d6 2 
dz 


(01-02)  + 


02 


_ 

Ps.aO  -  <t>2)h^z 


(13) 


where  Q_2  =  -Q.12  +  Qw2  +  O23  +  02ohm.  O23  =  /Cs.a ( 1  -  <P2)As 
(02  -  03)/[(?2+?3)/2],  Pi  =  PiRfToo/poo,  Ts,a  =  cPif/cv,s,a,  and  the 

subscript  s,a  and  s,c  mean  the  solid  part  of  the  anode  and  cathode, 
respectively. 

The  mass  balance  for  CV6,  yields  iho2out  —  hio2in  —  dio2  ar|d 
fiH2o  =  hn2o,out  =  hH2o,m  —  ho2.  The  mass  and  energy  balances 
for  CV6  state  that 


^  =  [o7  +  Hh2o(06)  -  Hh2o(07)  +  ^ox^y-(0ox  -  07)| 

Pox  Tox07,O 
Rf  Poxnc?c?7?z 

(10) 

where  Q7  =  -Q.67  +  Qw7  +  Q.7ohm.  Q67  =  h7As(  1  -  tpG)(67  -  d6), 
and  the  dimensionless  enthalpy  of  formation  of  a  compound  is 
defined  as  Hj  =  hiHi/(rhrefCpfroo)  where  the  subscript  i  refers  to 
a  substance  or  a  control  volume. 

2.4.2.  Anode  and  cathode  diffusion  layer  (CV2  and  CV6) 

The  electrode  is  divided  in  two  parts:  diffusive  layer  and  reactive 
layer.  The  purpose  of  the  diffusion  layer  is  to  provide  gas  uniformly 
to  the  reactive  layer  and  to  remove  the  byproducts  of  the  reaction. 
The  diffusion  layer  is  made  hydrophobic  to  avoid  the  electrolyte 
solution  to  go  to  the  gas  stream.  Therefore,  it  is  considered  that 
there  is  only  hydrogen  in  the  anode  diffusion  layer  (CV2)  and 
oxygen  and  water  (product  of  the  oxygen  reduction  reaction)  in  the 
cathode  diffusion  layer  (CV6). 

Diffusion  is  assumed  to  be  the  dominant  mass  transfer 
mechanism  in  the  diffusion  and  catalytic  layers.  As  the  friction  of 
the  substances  with  the  pore  walls  is  considered,  the  Knudsen 
flow  is  used  to  represent  the  fuel  and  oxidant  mass  flow  [23],  as 
follows: 

ii  =  -P(Pout-Pin)/%  =  2-  6  (11) 

where  D  =  B  [r  (8RT/7rM)1/,2  <pti]  is  the  Knudsen  diffusion  coeffi¬ 
cient,  r  the  pore  radius,  p  the  density,  R  the  universal  gas  constant, 


0.6  +  (07  -  06)  +  W(0s)h2o  -  H(06)h2o 

CP-f  J  (14) 

7s, c 
'  Ps.cC  ^ 

where  Q6  =  -Q56  +  QW6  +  Q67  +  Qeohm.  Ose  =  fes,c(l  -  <P6 )A 

(0s  -  0e)/[(?5  +  fe)/2],  and  ys>c  =  cPif/cv, s,c. 

2.4.3.  Anode  and  cathode  reaction  layer  (CV3  and  CV5) 

The  anode  (CV3)  is  assumed  to  contain  the  porous  electrode, 
fuel  (hydrogen)  and  electrolyte  aqueous  solution  (KOH)  that  is 
present  in  the  pores  of  the  electrode.  The  electrons  released  by  the 
oxidation  of  the  hydrogen  in  CV3  pass  through  an  external  circuit 
and  participate  in  the  reduction  of  the  oxygen  in  CV5.  The  cathode 
(CV5)  is  assumed  to  contain  the  porous  electrode,  oxygen  and  the 
electrolyte. 

The  Hydrogen  Oxidation  Reaction  (HOR)  and  the  Oxygen 
Reduction  Reaction  (ORR)  are  written  as  follows 

H2(g)  +  20H^aq)  -*2H2Oa)  +  2e  (15) 

lo2(g)  +  2H20(1)  +  2e-  -20H(aq)  +  H20(1)  (16) 

In  CV3  and  CV5,  the  mass  of  liquid  contained  in  these  CVs  is 
assumed  to  be  negligible  compared  to  the  solid  mass  of  the  elec¬ 
trode  for  the  thermal  analysis.  CV3  and  CV5  interact  thermally  with 
the  adjacent  control  volumes  by  conduction  and  the  with  the 
ambient  by  convection.  Heat  is  generated  by  the  electrochemical 
reactions  and  by  the  irreversible  mechanisms  (ohmic  resistance, 
activation  and  concentration  overpotential  losses). 
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Among  the  model  enhancements  of  the  present  study  are  the 
inclusions  of  the  heat  generation  due  to  activation  overpotential 
losses  (?ja)  and  concentration  overpotential  losses  (ijd.a).  However, 
the  ohmic  heat  generated  in  CV3,  CV4  and  CV5,  due  to  the  ion 
transport  resistance  of  the  solution  and  the  electric  resistance  of 
their  solid  components  has  been  previously  considered  for  the  AFC 
model  [1],  The  calculations  of  overpotential  losses  (r/a.  Vd,a  and 
?7ohm)  are  presented  in  the  electrochemical  model  section  of  this 
study. 

The  mass  and  energy  balances  for  CV3,  together  with 
the  anode  reaction  equation,  deliver  the  relations 
nH2  =  ™f/Mn2,  hH2o  =  2hH2,  mH2o  =  2nH2MH2o,  h0H~  =  2%,, 
hloH-  =  2hH2MOH-  and 

§=[a3-AH3+Ae3]  ^  (,7, 

where  Q3,act  =  ^aVi^eftpr^),  Q3dif  =  ^d,aty(*refCp,fr‘»)» 
0.3  —  -023  +  Ow3  +  034  +  O3ohm  +  Ch.act  +  03,dif.  and  Q.34  —  - 

(1  -  ¥>3) (03  -  04)2Asks,akmemb/(?4ks,a  +  ?3kmemb). 

Analogously,  the  mass  and  energy  balances  in  CV5,  and 
the  cathode  reaction  equation  deliver  hH20iin  =  2hH2o,out 
=  fio2,fi0H-ut  =  2ii02,  and 

^5  =  |q5  -  A Hs  +  AG51  -  ...  Ts’c.,  ,  e  (18) 

dr  L  Ps,c(l  -  ¥>5)?5?y?Z 

where  Qs.act  =  Vcl/(mrefCpjT„),  Q5,dif  =  Vd,CI/(mrefCp,fT«,), 
O5  =  -Q45  +  Qw5  +  O56  +  Osohm  +  05,act  +  Q5,dif.  and  CL45  =  - 
(1  -  <505) (^4  -  ^5)2AskSickmemb/(?4ks,c  +  fsftnemb)- 

Eqs.  (17)  and  (18)  account  for  the  reactions  dimensionless 
enthalpy  and  free  Gibbs  energy  variations  (AH3,AG3) 
=  ”h2 ( AH3 , AG3 ) / (thref cp f T 00 )  and  (AH5,AG5)  =  nH2(AH5,  AG5)/ 

(iflrefCpfToo). 

The  molar  enthalpies  of  formation  are  obtained  from  tabulated 
values  [26,27,28]  at  1  atm,  since  AH  is  independent  of  pressure,  and 
in  the  following  temperatures:  T2  for  H2(g),  T4  for  OHf-q)  and  T3  for 
H20(i).  For  CV5,  values  are  used  also  at  1  atm  and  at  Tg  for  02(g),  T4 
for  the  reactant  H20(i),  T5  for  the  product  H20(i)  and  OH^. 

The  molar  enthalpy  change  due  to  the  anode  reaction  is  given  by 

AH3  =  Eproducts[HH,(T,-)]  -  EreactantslHHKT,)]  and  We3  =  -AG3, 

AH3  is  the  CV3  reaction  enthalpy  change  (kj  kmol  1H2):  Vj  are  the 
stoichiometric  coefficients:  H,(T,)  is  the  molar  enthalpy  (kj  kmol-1) 
of  formation  at  a  temperature  Ti  of  reactants  and  products  of 
compound  i;  AG3  is  the  CV3  reaction  Gibbs  free  energy  change 
(kj  kmol-1H2),  and  We 3  is  the  maximum  (reversible)  electrical  work 
generated  due  to  the  reaction  in  CV3  (kj  kmol-1H2). 

The  reaction  Gibbs  free  energy  change  AG  is  a  function  of 
temperature,  pressure  and  concentrations, 

AG  =  AG°  +  RTlnQ  (19) 


The  molar  concentration  of  KOH,  [KOH],  is  obtained  from  the 
weight  percentage  of  KOH  in  solution,  y. 

For  the  cathode  (CV5),  AH  and  AG  are  obtained  similarly  to  the 
CV3  analysis,  obtaining  A H5  and  Wes  =  -AG5.  The  CV5  reaction 
quotient  is  Q5  =  [OH-]2/p^2,  in  which  p02  =  P6,out- 

2.4.4.  Electrolyte  and  membrane  (CV4) 

Typical  AFC  electrolytes  are  aqueous  solutions  of  KOH.  In  the 
AMFC,  the  electrolyte  comprises  an  aqueous  solution  of  KOH  sup¬ 
ported  by  an  inert  solid  medium.  In  the  prototype  shown  in  Fig.  1, 
a  cellulose-based  membrane  is  utilized. 

Eqs.  (15)  and  (16)  and  the  conservation  of  mass  in  CV4  require 
that  2hH2  =  hom, out  =  n0tr,in  =  2ho2,  2hH2  —  nH20in  =  fiH2o,out 
=  2ho2 .  In  CV4  the  mass  of  fluid  is  not  negligible  in  presence  of  the 
cellulose-based  membrane  mass  in  the  thermal  analysis.  Therefore, 
the  mass  and  energy  balances  for  CV4  state  that 

=  [0-4  +  Hoh  (85)  -  Hoh  (84)  +  Hh2o(03)  -  Wh2o(04)] 

Tmemb 

Pmemb?4?y?z 

(21) 

where  Q4  =  -Q34  +  Q.w4  +  Q45  +  CLtohm.  Pmemb  =  Psol9>4+ 

Pcel(l  -  <P4).  Cp  memb  =  Cp.sol^  +  Cp,cel(l  “  <P4).  and  kmemb 
=  ksol<P4  +  fccel(l  -  94)- 


2.5.  Electrochemical  model 

Based  on  the  electrical  conductivities  and  geometry  of 
each  compartment,  the  electrical  and  ionic  resistances,  (3(i2),  are 
given  by: 


i  =  1, 2, 6, 7{<p\ ,  <?7  =  0) 


(22) 


ft  =  - ,  i  =  3,4  and  5  (23) 

AsVj/3<Tj$j 

For  i  =  3,  4,  and  5,<r,-  =  <tsoi.  The  conductivities  of  the  diffusive 
layers,  02  and  <rg,  are  the  carbon-phase  conductivities  [29],  Finally, 
the  conductivities  of  CV1  and  CV7,  <j\  andff7,  are  given  by  the 
electrical  conductivity  of  the  bipolar  plate  material. 

The  dimensionless  potential  and  the  dimensionless  over¬ 
potential  are  defined  in  terms  of  a  given  reference  voltage,  Vref,  so 
that  Vt  =  Vi/Vre{  and  ft  _=  m/Vref. 

The  actual  potential  (V^)  provided  by  a  fuel  cell  results  from  the 
combination  of  anode  irreversible  potential  (Via),  cathode  irre¬ 
versible  potential  (VjjC)  and  the  ohmic  loss  (ifthm)  that  occurs 
throughout  the  fuel  cell  between  CV1  and  CV7,  as  follows: 


where  AG°  =  A H°  -  TAS°  is  the  standard  Gibbs  free  energy 
(kj  kmol-1H2);  A H°  is  the  standard  enthalpy  change  (kj  kmol-1); 
AS°  is  the  standard  entropy  change  (kj  kmol-1  K)  [gases  at  1  atm, 
25  °C,  species  in  solution  at  1  M,  where  M  is  the  molarity  =  (moles 
solute)/(liters  solution)],  and  Q  is  the  reaction  quotient.  In  the 
present  reaction  [Eq.  (15)]  the  resulting  expression  for  Q3  is 
Q3  =  { [OH(aq)]2pH2 }-  ' ,  where  [OH^q)]  is  the  molar  concentration 
of  the  alkaline  solution  and  ph2  =  P2,out.  i.e.,  the  H2  partial  pressure 
at  the  CV2  outlet. 

The  dissolution  of  potassium  hydroxide  in  water  is  represented 
by  the  following  chemical  equation: 

(20) 


V'i  =  V'i.a  +  V'ix  -  VoUm  (24) 

5ohm=^Eft  (25) 

The  reversible  electrical  potential  at  the  anode  and  cathode  are 
respectively  given  by  the  Nernst  equation,  as  follows: 

Vej  =  V° j  —  5^1nQj  (i  =  3,5  and  j  =  a,  c,  respectively) 

(26) 

where  Vej  =  AGi/(-nF),  Vej.  =  AG./(-nF). 


KOH(S)  — ►K(laq)  +  OH(aq) 
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At  the  anode  there  are  two  mechanisms  for  potential  losses:  (i) 
charge  transfer,  and  (ii)  mass  diffusion.  The  potential  losses  at  the 
anode  ( i)a )  and  cathode  (t]c)  due  to  charge  transfer  are  obtained 
implicitly  from  the  Butler-Volmer  equation  for  a  given  current  / 
[1,30,31],  as  follows: 


(i  =  3,5  and  j  =  a,c,  respectively) 
(27) 


where  a.j  are  the  anode  and  cathode  charge  transfer  coefficients, 
and  ioj  the  anode  and  cathode  exchange  current  densities, 
respectively. 

The  potential  losses  due  to  mass  diffusion  is  [1,30]: 


^=Jln  r- 


A',  wetium  j  J 


(i  =  3,5  and  j  =  a,c,  respectively) 
(28) 


The  limiting  current  densities  in  the  anode  and  cathode  (iumj) 
represent  borderline  cases  for  mass  transfer,  when  the  concentra¬ 
tion  of  reactants  in  the  catalyst  layer  located  at  the  interface  with 
the  gas  diffusion  layer  drop  to  zero  (P2,out  =  0  and  P6,out  =  0).  These 
are  the  maximum  currents  that  a  fuel  cell  can  provide.  For  higher 
current  values  the  system  collapses.  Then,  Eq.  (12)  is  re-written  as 
follows: 


iumj  =  =  2’ 6  and  J  =  a' c’  resPectively) 

(29) 

The  resulting  electrical  potential  at  the  anode  and  cathode  are 
respectively  given  by 

Vij  =  Vej  -  Vj  -  \vdj\  (j  =  a,c)  (30) 

The  absolute  values  of  r)dj-  are  used,  because  there  could  be 
Vdj  <  0  (cathodic  overpotential). 

2.6.  Fuel  cell  net  power  output 

In  order  to  estimate  the  available  power  provided  by  a  fuel  cell 
(IVnet),  it  is  necessary  to  subtract  the  energy  needed  to  supply  the 
fuel  cell  with  fuel  and  oxidant  (VVP).  Therefore  the  total  net  power 
(available  for  utilization)  of  the  fuel  cell  is 

Wnet  =  W-WP  (31) 

where  W  =  Vj  is  the  total  fuel  cell  electrical  power  output.  The 
dimensionless  pumping  power,  Wp,  is  given  by 

Wp  =  ^fSfJ'AP,  +  ^0XS0x^AP7  (32) 

where  S;  =  mrefr«,  Ri/VrefIref.  i  =  f,  ox. 

3.  Numerical  method 


Eqs.  (4),  (10),  (13),  (14),  (17),  (18),  (21),  for  dimensionless 
temperatures  and  Eq.  (12)  for  dimensionless  pressures,  and  the 
specified  initial  conditions  form  a  system  of  7  ordinary  differential 
equations  and  2  algebraic  equations.  The  unknowns  are  0j  and  Pi, 
i.e.,  the  temperatures  in  the  seven  control  volumes,  and  the  gas 
pressures  at  CV2  and  CV6  outlets.  Once  the  temperatures  and 
pressures  are  known,  the  electrical  potentials,  and  single  AMFC 
electrical  and  net  power  are  calculated  for  any  current  level. 


Two  numerical  methods  were  used.  The  first  method  calculates 
the  transient  behavior  of  the  system,  starting  from  a  set  of  initial 
conditions,  then  the  solution  is  marched  in  time  (and  checked  for 
accuracy)  until  a  steady  state  is  achieved  at  any  current  level.  The 
equations  are  integrated  in  time  explicitly  using  an  adaptive  time 
step,  4th-5th  order  Runge-Kutta  method  [32],  The  time  step  is 
adjusted  automatically  according  to  the  local  truncation  error,  which 
is  kept  below  a  specified  tolerance  of  10-6.  The  second  method  is  for 
the  steady  state  solution.  The  time  derivatives  are  dropped  from  Eqs. 
(4),  (10),  (13),  (14),  (17),  (18),  and  (21),.  Pressures  are  related  to 
temperatures  via  Eq.  (12).  The  system  reduces  to  seven  nonlinear 
algebraic  equations,  in  which  the  unknowns  are  the  temperatures  of 
the  seven  control  volumes.  This  system  is  solved  using  a  quasi- 
Newton  method  [32],  Convergence  was  achieved  when  the 
Euclidean  norm  of  the  residual  of  the  system  was  less  than  10  6. 


4.  Experiments 

4.1.  Experimental  set  up 

A  purpose-built  test  facility  was  implemented  in  the  laboratory 
for  conducting  fuel  cell  experiments,  as  shown  in  Fig.  la.  It  consists 
of  a  worktop  with  H2  and  O2  feeding  systems  with  respective 
manometers  and  flowmeters  mounted  on  the  wall,  an  adjustable 
low  power  electric  load,  built  with  nickel-chrome  alloy  wire,  and 
a  supporting  structure  for  single  fuel  cell  testing. 

A  prototype  of  a  new  single  AMFC,  which  includes  a  cellulose- 
based  membrane  impregnated  with  a  solution  of  water  and 
potassium  hydroxide  (KOH),  was  also  built-in  the  laboratory,  as 
shown  in  Fig.  lb.  The  membrane  is  composed  of  filter  paper  with 
a  specially  designed  composition  for  proper  electrolyte  absorption 
and  durability  [12],  providing  a  novel  and  inexpensive  method  of 
introducing  the  potassium  hydroxide  electrolyte  into  the  cell.  As 
shown  in  the  exploded  view  of  Fig.  3,  in  the  design  of  the  bipolar 
plates,  a  parallel  channel  configuration  was  utilized,  so  that  proper 
uniform  fuel  and  oxidant  gas  diffusion  is  achieved  starting  from  the 
backing  layer  of  the  electrodes.  The  membrane  is  sandwiched 
between  the  two  electrode  sheets  (anode  and  cathode).  The  solid 
electrolyte  or  the  cellulose-based  membrane  is  what  differentiates 
this  AMFC  from  other  AMFC  types.  Feedback  controlled  humidifiers 
containing  aqueous  potassium  hydroxide  solution  ensured  that  the 
membrane  was  not  dried  out  by  the  flow  and  kept  the  mass  fraction 
of  KOFI  in  the  electrolyte  at  a  constant  value  for  continuous 
operation. 


Fuel  Channels 


Fig.  3.  ExPloded  view  of  the  comPonents  of  the  single  AMFC  prototyPe. 
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The  AMFC  works  with  electrodes  composed  of  nickel  and  cobalt 
that  are  electrodeposited  onto  each  of  the  permeable  carbon  film 
sheets,  that  are  made  hydrophilic  on  the  reactive  side  that  face  the 
electrolyte  and  hydrophobic  on  the  side  facing  the  gas  channels. 
However,  in  order  to  perform  the  experimental  validation  of  the 
mathematical  model  for  this  study  as  a  base  case,  platinum-carbon 
electrode  sheets  were  used,  type  LT-250-EW,  from  BASF,  with  the 
load  of  5  g  m-2  of  30  wt.  %  Pt  supported  in  Vulcan  XC-72. 

4.2.  Experimental  procedure 

The  experimental  set  up  described  in  the  previous  section  was 
used  to  gather  the  experimental  data  set  required  for  the  experi¬ 
mental  validation  of  the  mathematical  model  introduced  in  this 
study.  Each  test  consisted  of  measuring  the  values  of  current  and 
voltage  in  order  to  quantitatively  describe  the  performance  of  the 
AMFC  by  polarization  and  power  curves.  This  was  done  by  varying 
the  load  resistance  shown  in  Fig.  la  from  open  circuit  (1  =  0)  to 
short  circuit  (maximum  I)  conditions. 

The  experimental  work  involved  the  acquisition  of  the  AMFC 
prototype  voltage  and  current  output  data  in  real  time.  This  task 
was  performed  through  the  utilization  of  a  computational  data 
acquisition  system  which  consisted  of  a  digital  multimeter  board, 
NI  PCI-4060,  a  NI  PCI-6703  analog  output  board,  and  a  SCXI-1127 
32-channel  high  voltage  multiplexer,  all  manufactured  by 
National  Instruments,  USA,  which  allowed  for  the  sequential  data 
acquisition  from  32  channels  at  interval  times  of  0.1  s.  All  the  data 
were  processed  by  a  suitable  software  application  to  convert  the 
signals  into  voltage,  current  and  power  readings. 

Measurements  of  current  and  voltage  were  taken  at  three 
different  times  of  the  experiment:  0, 15  and  30  minutes.  Using  the 
adjustable  load,  and  varying  the  load  resistance  progressively  so 
that  the  entire  cell  current  operating  range  was  covered  (0  <  /  < 
fmax).  it  was  possible  to  obtain  polarization  and  power  curves  for 
the  AMFC  prototype  in  the  three  time  settings.  Therefore  three  runs 
were  considered  for  each  membrane  KOH  concentration  corre¬ 
sponding  to  the  three  time  settings.  Additionally,  the  whole 
experimental  set  up  shown  in  Fig.  1  was  inside  a  temperature 
controlled  room,  in  which  the  air  temperature  was  set  to  25  °C. 
Steady  state  conditions  were  reached  after  30  min  before  starting 
all  the  experiments.  The  voltage  and  current  precision  limits  for 
each  experiment  were  computed  as  two  times  the  standard  devi¬ 
ation  of  the  3  runs  [33],  so  that  it  is  assumed  that  the  system 
response  follows  a  symmetric  unimodal  normal  distribution. 
Therefore,  it  is  admissible  to  state  that  95  %  of  the  system  responses 
exhibit  actual  voltage  and  current  equal  to  the  calculated  mean 
(from  the  sample  measurements)  plus  or  minus  twice  the  experi¬ 
mental  measurements  standard  deviation  [34],  Both  the  voltage 
and  current  bias  limits  were  negligible  in  presence  of  the  respective 
precision  limits.  Due  to  the  low  current  range  observed  in  the  AMFC 
prototype  (0  <  /  <  / max) ,  pumping  power  was  negligible  in  presence 
of  electrical  power,  therefore  VV„et  =  VV  =  VI.  As  a  result,  the 
measurements  uncertainties  to  produce  the  error  bars  in  the 
graphical  results  were  calculated  by: 


5.  Results  and  discussion 

In  this  section,  the  AMFC  prototype  behavior  is  studied  experi¬ 
mentally  and  numerically.  First,  the  cellulose-based  membrane 
electrolyte  solution  KOH  concentration  effect  on  AMFC  perfor¬ 
mance  is  analyzed.  Then,  the  experimental  measurements  are  used 
to  experimentally  validate  the  numerical  results  obtained  with  the 
mathematical  model  introduced  in  section  2,  using  5  (five)  different 
sets  of  measured  data  obtained  for  the  cellulose-based  AMFC 
shown  in  Fig.  1  for  several  electrolyte  concentrations.  The  proce¬ 
dure  consists  of  solving  the  inverse  problem  of  parameters  esti¬ 
mation  to  elaborate  an  empirical  correlation  for  the  electrodes 
exchange  current  density  as  a  function  of  electrolyte  concentration 
using  3  (three)  data  sets,  therefore  adjusting  the  mathematical 
model.  After  that,  the  conditions  corresponding  to  the  other  2  (two) 
data  sets  are  computationally  simulated  with  the  adjusted  model  to 
verify  the  agreement  between  numerical  and  experimental  results. 
Finally,  the  experimentally  validated  model  is  used  to  investigate 
the  AMFC  prototype  response,  temperature  distribution  and  pres¬ 
sure  drops. 

5.1.  The  effect  of  KOH  concentration  on  performance 

The  AMFC  prototype  of  Fig.  1  has  a  fuel  cell  square  section,  i.e., 
?y/?x  =  ?z/?x  =  3.95,  according  to  Fig.  2.  The  internal  configura¬ 
tion  is  given  by:;, /?x  =  ?7/?x  =  0.0789,  ?2/?x  =  ?e/?x  =  0.0108, 
?3/?x  =  fs/f*  =  0.0058,  ?4/?x  «  0.0197,  ?8/?x  =  ?9/?x 

=  0.3946,  ?t/?y  =  0.02,  and  ?c/?y  =  0.01.  The  bipolar  plates 
(interconnects)  comprise  two  parts:  the  gas  channels  (CV1  and 
CV7)  and  back  plates  that  also  act  as  mechanical  support,  which  are 
represented  by  L8  and  L9  in  Fig.  2. 

The  AMFC  response  operating  with  six  different  conditions, 
varying  the  KOH  mass  fraction  in  the  electrolyte  solution  (10  wt.  % , 
20  wt.  %  ,  30  wt.  %  ,  40  wt.  %  and  50  wt.  %  ),  is  shown  in  Fig.  4.  The 
measurements  of  current  and  voltage  were  taken  at  three  different 
times  of  the  experiment:  0,  15  and  30  minutes.  Since  the  AMFC 
prototype  performance  did  not  vary  significantly  on  these  three 
times  for  each  experiment,  it  is  reasonable  to  assume  that  the 
prototype  startup  transient  is  very  short,  which  is  probably  due  to 
the  relatively  small  size  of  the  cell,  thus  with  low  thermal  inertia. 

All  experimental  results  for  t  =  0  min  were  grouped  on  Fig.  5  and 
summarized  in  Table  1,  so  that  the  voltage  losses  that  occurred  and 
the  maximum  current  and  power  obtained  in  each  test  can  be 
compared.  The  maximum  power  was  obtained  for  yopt  =  40  wt.  % , 
which  shows  that  even  without  presenting  the  highest  initial 
potential,  the  cell  reached  the  highest  short  circuit  current  among 
all  tested  electrolyte  solution  mass  fractions,  and  showed  smaller 
potential  losses  than  with  the  other  concentrations,  which 
produced  the  highest  power  among  the  tested  KOH  mass  fractions. 
From  the  physical  point  of  view,  such  optimum  KOH  mass  fraction 
value  is  expected  based  on  the  analysis  of  two  extremes,  i.e.,  when 
y->0%,  hydroxyl  ions  become  unavailable  in  the  membrane, 
therefore  W  — >  0,  and  when  y  ->  1 00%,  the  flow  of  ions  decreases  due 
to  the  reduction  of  liquid  water  in  the  control  volume,  and  W  — ►  0  as 
well.  Therefore,  there  must  exist  an  intermediate  optimal  value  yopt, 
so  that  W  is  maximized. 

Fig.  6  summarizes  the  optimization  results  found  through  the 
experimental  measurements.  The  experimental  voltage  and  power 
measurements  corresponding  to  different  current  levels  are  plotted 
in  dimensionless  form  with  error  bars  calculated  according  to  Eqs. 
(33)  and  (34).  Both  maximum  voltage  and  power  were  found  for 
the  AMFC  prototype  operating  with  y  between  30  and  40  wt.  %  for 
three  current  levels  (I  —  0.56,  1.05,  and  2.52).  However,  as 
current  increases,  maximum  power  is  obtained  when  y  -» 40  wt.%. 
In  fact,  with  y  =  40  wt.%,  the  AMFC  reached  the  maximum  power 
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Fig.  4.  AMFC  prototype  experimental  potential  and  power  responses  for  five  different  electrolyte  solution  mass  fractions:  y  =  10  wt.  %  (a),  20  wt  X  (b),  30  wt.  %  (c),  40  wt.  %  (d)  and 
50  wt.  %  (e). 


1(A) 


Fig.  5.  Polarization  and  power  curves  experimentally  obtained  with  different  elec¬ 
trolyte  mass  fraction  (y)  values  at  t  =  0  min. 


and  operating  current  in  short  circuit  among  all  tested  cases, 
indicating  that  the  maximum  performance  is  obtained  for 
yopt  =  40  wt.%  according  to  the  experimental  measurements. 

The  effect  of  KOH  concentration  on  alkaline  fuel  cells  has  been 
investigated  theoretically  before  by  Zhang  et  al.  [35]  that  conducted 
a  multi-objective  optimization  of  a  typical  alkaline  fuel  cell  using 


Table  1 

Experimental  values  for  open  circuit  and  short  circuit  at  t  =  0  min. 


y  (wt.  %  KOH)  i  (A)  V(V) 

"To  000  0.92 

2.56  0.54 

20  0.00  0.94 

3.28  0.62 

30  0.00  0.98 

3.48  0.73 

40  0.00  0.93 

4.11  0.72 

50  0.00  0.98 

3.63  0.63 


W(W) 

0.00 

1.39 

0.00 

2.05 

0.00 

2.55 

0.00 

2.98 

0.00 
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a  steady  state  alkaline  fuel  cell  model.  The  authors  reported  an 
algebraic  correlation  to  obtain  the  optimal  KOH  molar  concentra¬ 
tion  for  maximum  alkaline  fuel  cell  performance  as  a  function  of 
operating  temperature.  The  present  study  corroborates  experi¬ 
mentally  the  theoretically  expected  trend,  although  for  a  different 
alkaline  fuel  cell  type,  i.e.,  the  herein  proposed  cellulose-based 
membrane  alkaline  fuel  cell.  Therefore,  the  optimal  KOH  concen¬ 
tration  for  maximum  power  output  is  a  fundamental  feature  that 
should  be  pursued  in  any  alkaline  fuel  cell  type. 

5.2.  Experimental  model  adjustment  and  validation 

A  computational  code  was  written  in  Fortran  language  based  on 
the  numerical  method  described  in  section  3  to  obtain  the  solution 
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Fig.  6.  Experimental  AMFC  potential  (a)  and  power  (b)  output  values  as  functions  of 
KOH  mass  fraction  for  three  current  levels  (1  =  0.56, 1.05  and  2.52). 


to  the  AMFC  mathematical  model,  i.e.,  temperatures  and  pressures, 
and  then  computing  AMFC  potentials  and  net  power  for  the  entire 
current  operating  range  (0  <  I  <  /max).  The  input  parameters  are  the 
geometric  features  and  physical  properties  of  the  single  AMFC.  The 
geometry  and  several  physical  properties  were  measured  directly 
from  the  AMFC  prototype,  and  are  listed  in  Table  2.  The  remaining 
necessary  data  were  taken  from  the  technical  literature 
[1,14,26,28,29,36-40],  The  reference  case  (geometry  and  proper¬ 
ties)  for  all  simulations  in  this  study  is  the  AMFC  prototype. 

The  current  (!)  is  considered  an  independent  variable  in  the 
numerical  simulation.  The  fuel  cell  electrical  and  net  power  were 
computed  by  starting  from  open  circuit  (7  =  0),  and  proceeding 
with  current  increments  (A/  =  0.1)  until  the  electrical  current 
reached  the  values  experimentally  obtained  with  AMFC  prototype. 

As  pointed  out  earlier  in  the  text,  the  model  experimental 
validation  was  performed  using  5  (five)  different  sets  of  measured 
data  obtained  with  different  KOH  mass  fraction  in  the  electrolyte 
solution  (10  wt.  % ,  20  wt.  % ,  30  wt.  % ,  40  wt.  %  and  50  wt.  %  ).  There 
are  some  parameters  that  change  with  solution  concentration 
(ksoi,  (0Soi  and  crS0|)  and  the  appropriate  values  were  taken  from  the 
literature  [38,39]  and  are  shown  in  Table  3.  For  clarity,  the  geometry 
of  the  AMFC  prototype  used  in  the  experiments  is  summarized  in 
Table  4.  The  model  experimental  validation  initial  step  was  to 
perform  the  model  adjustment. 

The  model  adjustment  procedure  consisted  of  solving  the 
inverse  problem  of  parameters  estimation  (IPPE)  [41],  using  the 
mathematical  model,  i.e.,  by  turning  a  variable  into  a  parameter, 
and  vice-versa.  For  that,  what  was  originally  a  variable  is  imposed 
to  the  model  as  an  input  parameter.  In  this  study,  the  variables 
available  from  the  experiments  were  the  measured  voltage  and 
current  output  of  the  AMFC  prototype.  Therefore,  two  parameters 
could  become  variables  in  the  model.  The  procedure  started  by 
selecting  the  parameters  to  become  the  variables  to  be  determined. 
The  anode  and  cathode  exchange  current  densities  were  selected 
due  to  their  expected  dependency  on  electrolyte  solution  concen¬ 
tration,  which  was  the  parameter  that  was  experimentally  varied. 
The  exchange  current  density  determines  the  magnitude  of  the 
oxidation  and  reduction  reactions  that  occur  in  the  electrodes  at 
the  equilibrium  potential.  Even  though  the  net  current  is  zero  at 
equilibrium,  balanced  faradaic  activity  happens,  and  is  expressed  in 
terms  of  exchange  current  density  [30,31]  as  follows: 

to  =  nFk°QCg  “>  (35) 

where  k°  is  the  standard  rate  constant,  and  Co  and  Cr  the  concen¬ 
tration  of  oxidized  and  reduced  states  respectively. 


B  =  0.156  [42] 

cp,f=  14.307  kj  kg1  K1 

Cp.ox  =  0-918  kj  kg-1  K  ' 

cv,f  =  10.183  kj  kg1  K_1 
cv  ox  =  0.658  kj  kg-1  K"1 

fret  =  1  A 

kf=  0.182  Wm1  K1 
kox  =  0.0267  W  m1  K  ' 
fcp  =  0.1298  W  m1  K-1 
ks,a  =  ks,  c  =  0.1  W  m-1  K_1 
ksoi  =  0.571  W  m-1  K-1 
K2,  K6  =  4.2  x  1014  m2  [14] 
K3>  K5  =  4.2  x  10-16  m2  [14] 
pf=  0.127  MPa 
Pox  =  0.134  MPa 
poo  =0.1  MPa 
q  =  2.5 


Vref=  1  V 
ac  =  0.5 

Cl,  Cv  =  2 

p,  =  8.96  x  10  6  Pa  s 

ft 7  =  20.7  x  10  6  Pa  s 
tr sol  =  53.2  O'1  m-1 
ff,,  (17  =  1.5  X  107  fi-1  nr1 
o2,  o6  =  8570  fi-1  m-1 
02.  06  =  0.0085 
03,  05  =  0.172 
04  =  0.71 


Table  2 

Physical  properties  used  as  reference  case  in  the  numerical  solution. 

Rf  =  4.157  kj  kg-1  K-1 
Rox  =  0.2598  kj  kg-1  K  1 
tref  =  10“3  s 
Tfl  Tox,  Too  =  298.15  K 
Lfwi  =  50  Wm“2  K1,  i  =  1  to  7 
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Table  3 

Concentration-dependent  parameters. 


y  10  wt.% 

40  wt.% 

50  wt.% 

fcsoi  0.61 

0.57 

0.60 

Pso\  1040.84 

1406.39 

1532.69 

ffsoi  20.08 

53.22 

47.61 

i0,a  0.0021 

0.0058 

0.0046 

i0,c  0.0002 

0.0019 

0.0005 

Table  4 

Geometry  of  the  AMFC  prototype  of  Fig.  1. 

VT  =  2.5  x  10-4m3 

6/fx 

=  fe/k  =  0.0058 

fy/fx  =  fz/fx  =  3.95 

W?x 

=  0.0197 

fr/fx  =  ?7/?x  =  0.0789 

h/h 

=  £9/£x  =  0.3946 

?2/?x  =  f6/f*  =  0.0108 

ft/fy 

=  0.02 

fc/fy 

=  0.01 

Eq.  (35)  clearly  shows  that  if  the  electrolyte  solution  concen¬ 
tration  varies,  new  equilibrium  conditions  will  be  established, 
therefore  a  new  exchange  current  density  will  result.  Based  on  that, 
it  is  herein  proposed  to  elaborate  empirical  correlations  for  the 
electrodes  exchange  current  densities  as  functions  of  electrolyte 


solution  concentration  using  3  (three)  data  sets,  corresponding  to 
y  =  10  wt.  %  ,  40  wt.  %  and  50  wt.  %  .  The  numerical  solution  of  the 
IPPE  was  obtained  for  the  steady  state  at  each  current  incremental 
value,  and  delivered  exchange  current  values  that  showed  a  strong 
dependency  on  current  for  0  <7  —  1,  which  stabilized  in  the  range 
1  ~J  <  fmax  for  the  3  tested  electrolyte  solution  concentrations.  This 
effect  indicates  that  it  is  reasonable  to  expect  that  the  exchange 
current  density  for  1  ~I  <  7max  could  be  empirically  correlated  as 
a  function  of  the  electrolyte  solution  concentration. 

As  a  result,  the  AMFC  prototype  response  was  then  numerically 
simulated  using  the  stabilized  values  of  exchange  current  densities 
shown  in  Table  3  for  y  =  10  wt.  %  ,  40  wt.  %  and  50  wt.  %  .  Fig.  7 
shows  in  solid  lines  the  numerically  calculated  potential  curves 
for  all  tested  electrolyte  solution  concentrations  in  the  range  0  <  I  < 
fmax,  and  the  experimentally  measured  voltage  and  net  power  with 
error  bars  that  were  calculated  with  Eqs.  (33)  and  (34).  It  is  noted 
that  the  experimental  uncertainties  both  for  potential  and  power 
data  were  small  for  low  currents,  and  increased  as  current 
increased.  For  1  ~f  <  fmax,  the  numerical  results  lie  within  the  error 
bars  as  expected,  therefore  with  good  qualitative  and  quantitative 
agreement  with  the  experimental  results.  Regarding  the  power 
curves,  there  was  good  agreement  both  quantitatively  and  quali¬ 
tatively  between  the  simulation  results  and  the  experimental  data. 
The  numerical  results  remained  within  the  error  bars  throughout 
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Fig.  7.  Experimental  model  adjustment  using  y  =  10  wt.  %  (a),  40 ' 


;  (b),  and  50  wt  %  (c). 


EM.  Sommer  et  al.  /  Journal  of  Power  Sources  213  (2012)  16—30 


the  tested  current  range  (0  <  /  <  4.5)  and  the  maximum  AMFC 
prototype  measured  power  was  W =3. 

The  exchange  current  density  defines  the  activity  of  the  elec¬ 
trode  and  strongly  influences  the  fuel  cell  performance.  Using  the  3 
(three)  (io,a,  io,c)  values  of  Table  3,  the  following  empirical  correla¬ 
tions  were  obtained  for  the  electrodes  used  in  the  AMFC  prototype 
shown  in  Fig.  1 : 

!'0,a  =  -2.14  x  10~7y3  +  1.53  x  l(T5y2  -  0.00019  y  +  0.00272 

(36) 

i0c  =  -4.35  x  10-9y4  +  3.11  x  l(T7y3  -5.16  x  l(T6y2 
+  3.08  x  10-5y  -  0.00014 


After  performing  the  model  adjustment,  the  next  step  was  the 
model  experimental  validation.  This  task  was  achieved  by  using  the 
remaining  2  (two)  experimental  data  sets,  i.e.,  the  measured 
voltage  and  current  obtained  with  y  =  20  wt.  %  and  30  wt.  %  .  Eqs. 
(36)  and  (37)  were  then  used  to  determine  (io,a,  io,c)  for  y  —  20  wt.  % 
and  30  wt.  %  .  The  calculated  (to,a>to,c)  values  were  then  used  as 
input  to  the  mathematical  model  in  order  to  obtain  numerically  the 
AMFC  prototype  polarization  and  power  curves.  The  numerically 
simulated  AMFC  response  is  shown  in  solid  lines  in  Fig.  8  fory  =  20 
wt.  %  and  30  wt.  %  .  The  voltage  curves  lie  within  the  error  bars  for 
1  ~J  <  /max,  and  the  power  curves  for  the  entire  tested  range. 
Therefore,  there  was  good  agreement  between  numerical  and 
experimental  results  with  the  use  of  the  exchange  current  density 
empirical  correlations,  so  that  the  model  could  be  reliably  used  to 
simulate  and  analyze  the  AMFC  response  under  different  geometric 

(37) 

and  operating  conditions  for  the  range  10  wt.%  <  y  <  50  wt.%. 


I 
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Fig.  8.  Experimental  model  validation  using  y  =  20  wt  %  (a),  and  30  wt  %  (b). 


5.3.  AMFC  simulation 

In  order  to  illustrate  the  application  of  the  experimentally 
validated  AMFC  mathematical  model,  numerical  simulations  are 
conducted  fory  =  40  wt.%,  since  it  was  shown  experimentally  that 
this  is  the  optimal  electrolyte  solution  for  maximum  AMFC  net 
power  output.  Fig.  9  shows  the  resulting  polarization  and  net 
power  curves  obtained  with  the  simulation.  The  actual  open  circuit 
potential  is  equal  to  the  reversible  potential,  because  in  the  math¬ 
ematical  model  it  is  assumed  that  there  are  no  losses  due  to  species 
crossover  between  the  two  electrodes  and  from  internal  currents. 
Therefore,  with  the  model,  it  is  possible  to  verify  the  contribution  of 
anode  and  cathode  potentials  separately  to  the  resulting  total  AMFC 
potential  Vi.  As  expected,  cathode  potential  is  smaller  than  anode 
potential. 

Fig.  10  shows  the  AMFC  internal  temperature  distribution.  The 
temperature  along  the  AMFC  increases  as  current  increases,  since 
more  heat  is  being  generated  by  electrochemical  reactions  and 
potential  losses  (ohmic  resistance  and  overpotentials  -  activation 
and  diffusion).  The  CV5  shows  a  slightly  higher  temperature  than 
CV3  because  it  is  the  control  volume  in  which  the  oxygen  reduction 
reaction  takes  place,  which  releases  more  heat  than  the  hydrogen 
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Fig.  9.  The  AMFC  prototype  polarization  and  power  curves  obtained  numerically  with 
the  experimentally  validated  model  for  y  =  40  wt  %  . 
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Fig.  10.  The  thermal  response  of  the  AMFC  prototype  obtained  numerically  with  the 
experimentally  validated  model  for  y  =  40  wt  %  . 
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Fig.  12.  The  dynamic  thermal  response  of  the  AMFC  prototype  according  to  changes  in 
current  load  requirement  obtained  numerically  with  the  experimentally  validated 
model  for  y  =  40  wt.  %  . 


oxidation  at  CV3.  These  simulation  results  show  that  temperature 
varies  considerably  as  fuel  and  oxidant  flow  across  the  cell  during 
operation,  therefore  the  commonly  used  uniform  temperature 
assumption  for  fuel  cell  operation  only  holds  for  low  current  values. 
Such  temperature  spatial  dependence  is  important  to  be  consid¬ 
ered  in  order  to  obtain  representative  AMFC  power  output  values 
since  exchange  current  densities  and  other  parameters  depend  on 
temperature.  Such  temperature  spatial  dependency  is  expected  to 
be  more  pronounced  when  a  larger  current  interval  is  considered  or 
for  AMFC  stack  design. 

The  assessment  of  pressure  drops  in  the  gas  channels  (CV1  and 
CV7)  allows  for  obtaining  the  actual  value  of  output  power  produced 


Fig.  11.  The  dimensionless  pressure  drops  in  the  gas  channels  of  the  AMFC  prototype 
obtained  numerically  with  the  experimentally  validated  model  for  y  =  40  wt.  %  . 


by  the  AMFC  via  subtraction  of  the  power  required  for  pumping  the 
fuel  and  the  oxidant  from  the  electrical  power,  as  stated  by  Eq.  (31 ). 
The  simulation  results  of  Fig.  11  show  that  CV7  pressure  drop  is 
higher  than  in  CV1,  since  oxygen  is  roughly  ten  times  more  viscous 
than  hydrogen.  For  low  currents,  pressure  drops  are  small,  since 
small  fuel  and  oxidant  mass  flow  rates  are  needed.  However,  as 
current  increases,  pumping  power  should  not  be  neglected. 

Finally,  as  an  example  of  the  model  application  in  a  transient 
simulation,  a  transient  analysis  is  conducted  numerically  with  the 
model.  Fig.  12  shows  the  thermal  transient  response  of  the  AMFC 
when  the  load  current  requirement  changes  from  /  =  4  to  J  =  2, 
then  back  to  /  =  4.  Besides  information  on  the  time  required  to 
achieve  new  current  levels,  the  results  demonstrate  that  temper¬ 
ature  varies  significantly  from  one  current  level  to  another,  an  effect 
that  becomes  increasingly  important  as  current  increases,  due  to 
the  increase  in  fuel  cell  heat  generation  at  the  reacting  sites.  Note 
also  the  temperature  gradients  within  the  cell,  i.e.,  between  the  gas 
channels  and  membrane  electrodes  assembly.  Such  fuel  cell  oper¬ 
ating  temperature  variation  is  commonly  disregarded  in  most  fuel 
cell  studies  found  in  the  literature,  normally  under  the  assumption 
that  the  cooling  system  is  capable  of  equilibrating  fuel  cell 
temperature  at  any  desired  setpoint  independently  of  current, 
which  is  true  only  for  the  cooling  channels,  but  not  for  the  internal 
fuel  cell  components,  as  the  present  model  demonstrates. 


6.  Conclusions 

In  this  paper,  an  AMFC  dynamic  mathematical  model  with 
temperature  dependence  on  space  and  current  was  introduced  and 
experimentally  validated  by  direct  comparison  with  output  voltage 
and  power  measurements  performed  in  a  cellulose-based  AMFC 
prototype  in  the  laboratory.  The  single  AMFC  voltage  and  power 
output  were  obtained  numerically  with  the  mathematical  model, 
and  the  results  show  good  qualitative  and  quantitative  agreement 
with  the  measured  experimental  data.  Appropriate  dimensionless 
groups  were  identified  to  formulate  the  model,  and  the  results 
presented  in  normalized  charts  for  general  application. 
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It  is  concluded  from  the  experimental  results  that  the  startup 
transient  is  short  and  that  there  are  optimal  values  of  KOH  mass 
fraction  y  (~40  wt.  %  )  which  lead  to  maximum  power,  that  are 
herein  shown  experimentally  for  the  first  time.  The  model  was  used 
to  formulate  empirical  correlations  for  the  exchange  current 
density  (io)  in  the  electrodes  with  respect  to  the  electrolyte 
concentration,  and  could  be  used  for  future  AMFC  development. 

Temperature  gradients  along  the  fuel  and  oxidant  flow  path  are 
shown  to  exist  and  increasingly  more  important  as  operating 
current  increases,  both  in  steady  state  or  dynamic  operation, 
mainly  when  the  cell  is  subject  to  changes  in  operating  current 
regime  in  time,  which  is  a  frequent  demand  in  any  fuel  cell  appli¬ 
cation.  For  AMFC  stacks,  significant  temperature  spatial  gradients 
are  also  expected  since  large  coolant  flow  rates  are  limited  due  to 
induced  high  pressure  drops  in  the  cooling  channels,  and  should  be 
taken  into  account  in  optimal  AMFC  stack  design.  Therefore,  the 
obtained  results  demonstrate  the  importance  of  considering 
temperature  dependence  on  space  and  current  in  AMFC  mathe¬ 
matical  models.  The  model  presented  in  this  paper  accounts  for 
temperature  gradients  in  the  flow  direction  and  dependence  on 
operating  current.  The  model  also  considers  the  pressure  drops  in 
the  gas  channels,  so  that  AMFC  net  power  is  computed  for  more 
realistic  results.  As  a  result,  it  is  reasonable  to  state  that  the 
combination  of  accuracy  and  low  computational  time  allow  for  the 
future  utilization  of  the  model  as  a  reliable  tool  for  AMFC  simula¬ 
tion,  control,  design  and  optimization  purposes. 
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